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Abstract 

We investigate the quantum tunneling of a Bose-Einstein condensate con- 
fined in a optical trap. We show that periodic pulses of coherent matter are 
emitted from the trap by using an oscillating energy barrier. Moreover, the 
emitted fraction of condensed atoms strongly increases if the period of oscil- 
lation of the height of the energy barrier is in parametric resonance with the 
period of oscillation of the center of mass of the condensate inside the poten- 
tial well. Our model is analyzed by numerically solving the nonpolynomial 
Schrodinger equation (NPSE), an effective one-dimensional equation which 
describes the macroscopic wavefunction of Bose condensates under transverse 
harmonic confinement. The range of validity of NPSE is discussed and com- 
pared with that of Gross-Pitaevskii equation. 
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I. INTRODUCTION 



Nowadays dynamical properties and resonant phenomena of Bose-Einstein condensates 
are the subject of many experimental and theoretical invertigations. Recently, the para- 
metric excitation of cold trapped atoms in far-off-resonance optical lattices has been exper- 
imentally obtained by modulating the potential depth [1,2]. Another kind of parametric 
resonance, namely the parametric resonance between the collective oscillations of a trapped 
Bose-Einstein condensate and the oscillations of the confining harmonic trapping potential 
has been demonstrated with numerical simulations [3]. Moreover, by imposing particular 
values of anisotropy to the trapping harmonic potential, it has been shown that different 
collective modes of the condensate are in resonance and the oscillations can be chaotic [4,5]. 

In this paper we discuss a mechanism to generate coherent matter pulses by means of 
macroscopic quantum tunneling of a Bose condensate through the barrier of a potential well 
[6]. In our model the Bose-Einstein condensate is confined in the vertical axial direction 
by two Gaussian optical barriers and in the transverse direction by a magnetic or optical 
harmonic potential. The zero-temperature dynamics of Bose-Einstein condensate under 
transverse harmonic confinement is well described by the one-dimensional nonpolynomial 
Schrodinger equation (NPSE) we have recently derived from the three-dimensional Gross- 
Pitaevskii equation (3D GPE) [7,8]. In the first part of the paper we discuss the range 
of validity of the 3D GPE and that of the NPSE. In the second part of the paper by the 
numerical integration of NPSE we show that it is possible to generate periodic waves of 
coherent matter by periodically changing the height of the lower-lying Gaussian barrier. 
Moreover, the emission probability is strongly enhanced when the oscillation of the height 
of the energy barrier is in parametric resonance with the oscillation of the center of mass of 
the condensate inside the potential well. 
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II. CONDITIONS ON THE GROSS-PITAEVSKII EQUATION 



The time-dependent three-dimensional Gross-Pitaevskii equation (3D GPE) [9], that 
describes the macroscopic wavefunction (or order parameter) ip(r, t) of the Bose condensate, 
is given by 



d 

tk—ij(r,t) 



.^-V 2 + U(r)+gNmr,t)\ 2 



ri>(r,t), (1) 



where U(r) is the external trapping potential and g = AixTi 2 a s /m is the scattering amplitude 
and a s the s-wave scattering length. N is the number of condensed bosons and the wave- 
function is normalized to one. Note that the 3D GPE is a nonlinear Schrodinger (Hartree) 
equation and the nonlinear term is due to the interatomic interaction described by the 
two-body pseudo-potential V(r, r') = gS^(r — r') [10]. 

The 3D GPE is accurate to describe a condensate of dilute gas of bosons only near zero 
temperature, where thermal excitations can be neglected [10]. The dilute gas condition is 
given by a 3 p <C 1, where a s is the s-wave scattering length of the interatomic potential and 
p is the density of the Bose gas. This diluteness condition means that the range (a s ) of the 
interatomic potential is much smaller than the average distance (p -1 ^ 3 ) between atoms. By 
introducing the healing length 

£ = , (2) 

^8na s p 

the dilute gas condition can be written as a s <C £. Note that the condensate can be dilute 
but strongly interacting if its interatomic energy (4irh 2 a s p/m) is much larger than its kinetic 
energy. In the case of an axially symmetric harmonic confinement by using the healing length 
£ the dilute and strongly- interacting gas condition can be written as a s < ^ < a^, where 
an = \f^~ is the characteristic length of the harmonic confinement with ujh = (^x^z) 1 ^ 3 - 
Another interesting condition is the one-dimensional gas condition. One has the ID 
regime when the potential energy of the transverse confinement is much larger than the 
interatomic energy of the Bose condensate. In the case of a transverse harmonic confinement, 
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again using the healing length £, the ID gas condition for a Bose condensate can be written 



as a_i_ <C £, where aj_ = is the characteristic length of the transverse harmonic 

confinement. It follows that, if the inequalities a s C ajj < ai C ^ are satisfied, then the 
Bose condensate is dilute, cigar-shaped and one-dimensional. 

It is important to observe that for £ am, where am is the ID scattering length given by 
a \D = ci\/a s [10], the Bose gas enters in the Tonks- Girardeau regime: a one-dimensional and 
very dilute gas of impenetrable Bosons for which the Bose-Einstein condensation is absent 
[11]. Thus in the Tonks- Girardeau regime the 3D GPE cannot be applied. In conclusion, 
the 3D GPE is valid if the healing length £ of the Bose condensate satisfies the condition 

a s <£<ai/a s , (3) 

and it describes a one-dimensional Bose condensate if a s <C aj_ <C £ < a\/a s . 

III. NONPOLYNOMIAL SCHRODINGER EQUATION 

In many experiments with Bose-Einstein condensates the external trapping potential can 
be described by a harmonic potential in the transverse direction and a generic potential in 
the axial direction: 

U(r) = \nuj\{x 2 + y 2 ) + V(z) . (4) 

This external potential suggests to map the 3D GPE into an effective ID equation, which 
simplifies greatly the solution of the 3D GPE. This problem is not trivial due to the nonlin- 
earity of the GPE. 

The 3D GPE is the Euler-Lagrange equation of the following Lagrangian density 

d h 2 



iP(r,t). (5) 



L = r(r,t) 

For the wavefunction we choose the following variational ansatz 

V>(r, t) = 4>{x, y, t; a(z, t)) f(z, t) , (6) 



where both and / are normalized to one and is represented by a Gaussian: 

(f)(x,y,t;a(z,t)) 



g 2a(z,ty 
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Moreover we assume that the transverse wavefunction is slowly varying along the axial 
direction with respect to the transverse direction, i.e. V 2 ~ V^_0 where = ^ + Jp. 
By inserting the trial wave-function in the Lagrangian density and after spatial integration 
along x and y variables the Lagrangian density becomes 
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The Euler-Lagrange equations with respect to /* and a read 
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The second Euler-Lagrange equation reduces to an algebric relation providing a one to one 
correspondence between a and /: a 2 = a\\Jl + 2a s N\f\ 2 , where a± = \J * s ^ ne oscma tor 
length in the transverse direction. One sees that a depends implicitly on z and t because of 
the space and time dependence of \f\ 2 . Inserting this result in the first equation one finally 
obtains 

gN |/| 2 



dt J 



2m dz 2 
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This equation is a time-dependent non-polynomial Schrodinger equation (NPSE). 

We observe that under the condition a s iV|/| 2 <C 1 one has a 2 = a\ and NPSE reduces 

to 

h 2 d 2 



lU di f 



2m dz 2 



V + 



2vra 2 1/1 



(12) 



where the additive constant huj± has been omitted because it does not affect the dynamics. 
This equation is a fD Gross-Pitaevskii equation (ID GPE). The nonlinear coefficient g' 
of this ID GPE can be thus obtained from the nonlinear coefficient g of the 3D GPE by 
setting g' = g/(2Tra 2 L ) such that V(z, z') = g'S(z — z') = (2h 2 /maiD)5(z — z') is the effective 
ID interatomic pseudo-potential and the effective ID scattering length is air, = a±/ a s, in 
agreement with the definition of the previous section. Note that the limit a s N\f\ 2 <C 1 is 
precisely the condition for the one-dimensional regime where the healing length £ is larger 
than a±. But, as shown in the previous section, the 3D GPE and consequently the NPSE 
are no more valid if £ a\/a s . It follows that the ID GPE is valid only where the condition 
a 2 /a\ < a s N\f\ 2 1 is satisfied. 

In papers [9,10] we have tested the accuracy of the NPSE in the determination of the 
ground-state and collective oscillations of the condensate with axial harmonic confinement 
and also in the description of tunneling through a Gaussian barrier. In particular, we have 
compared NPSE with the the full 3D GPE. The conclusion of these investigations is that 
NPSE is very accurate in the description of cigar-shaped condensates (a z /a± > vTO) both 
in the 3D regime and in the ID regime (a s <C £ < a\/a s ). 

Note that the NPSE has been recently used by Massignan and Modugno to study the 
dynamics of a Bose condensate in a optical lattice [12]. 

IV. PERIODIC EMISSION OF MATTER WAVES 

In this section we use the NPSE to study a mechanism which produces periodic emission 
of matter waves by means of the quantum tunneling of a Bose condensate. We consider a 
cigar-shaped condensate confined by a harmonic potential in the transverse direction and 
under the action of the gravity potential mgz and two Gaussian functions that model a 
confining potential well in the vertical axial direction: 

V(z) = V 1 e -(^i) 2 /- 2 + V 2 e -(^) 2 /- 2 + mgz . (13) 
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Such a configuration can be experimentally obtained by using two blue-detuned laser beams 
(perpendicular to the axial direction) which are modelled by the two Gaussian potentials. 
By varying the intensity of the lower-lying laser beam one controls the height of the energy 
barrier and therefore the tunneling probability. The periodic height given by: 

2tv 

V 1 (t) = V a + V b sin(—t), (14) 

T 

where r is the period of oscillation. 

The NPSE is integrated by using a finite-difference predictor-corrector method [13,14,15]. 
The initial wave function of the condensate is found by solving the equation with imaginary 
time. In Figure 1 we plot the axial density profile p(z) of the matter-waves coming out from 
the potential well with r — 1.5. The figure shows that, as expected, the emission of pulses is 
periodic. Moreover, both the period of emission and the tunneling fraction Pt depend on r. 
In the inset of Figure 2 the tunneling probability Pt is shown as a function of r for a Bose 
condensate with N = 10 3 and iV = 10 4 atoms. The emission probability has its absolute 
maximum near r = 0.8, that is the average period of the oscillations of the center of mass 
of the condensate inside the potential well. 

The enhancement of the emission probability of coherent matter near r = 0.8 has a 
classical explanation. We have verified [6] that it corresponds to the condition of resonance 
between the oscillation of the energy barrier and the oscillation of the center of mass of 
the condensate. Close to the resonance condition the phenomenon of coherent emission of 
atoms is no more due to quantum tunneling but to the parametric resonance [16] between 
the period of oscillation of the center of mass of the condensate and the period of oscillation 
of the energy barrier. This might explain the different N dependence of Pt and P T as shown 
in Figure 2: Pj, is dominated by quantum effects and, as such, is enhanced by interaction, 
while Pt is essentially a classical effect and is almost N independent, in agreement with the 
numerical data shown in Figure 2. 

Note that in Figure 2 there are other local maxima whose positions follow the text-book 
[16] resonant condition formula r = (n/2) r osc , where r osc is the period of oscillation of the 
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center of mass in the absence of an oscillating barrier, with n an integer number. 



CONCLUSIONS 

In this paper we have studied the controlled emission of coherent matter pulses from a 
trap by changing the period of oscillation of the energy barrier. The experimental inves- 
tigation of the physical configuration we have considered in this paper may contribute to 
the realization of novel phenomena in atom lasers. Finally, we observe that the parametric 
driving of Bose-Einstein condensates can be obtained by current experiments using optical 
dipole forces with far-detuned laser beams. By varying the intensity of the laser beams one 
can control the height of the energy barrier that confines the condensed sample. 
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FIG. 1. Axial density profile of Bose condensed 23 Na atoms tunneling through the potential 
well (Eq. (13)) with oscillating barrier (Eq. (14)): V a = 300 and V& = 200. Results obtained by- 
solving NPSE. Period of oscillation: r = 1.5. Effective chemical potential of the initial condensate: 
fx eff = 57.43. Number of 23 Na atoms: N = 10 4 . Scattering length: a s = 30 A. Note that 
Pt = 0.26 with V\{t) = V a — 14 = 100. Length z in units a z = (h/muj z ), where u z = cjj_/10 with 
uj± = 2ir kHz; energy in units hto z and time in units UJ^ 1 - 
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FIG. 2. Tunneling ratio Pt/-?? at t = 5 as a function of the period r of oscillation of the 
lower-lying Gaussian barrier (Eq. (13,14)), as obtained by solving NPSE, with V a = 300 and 
V b = 200 (Eq. (6)). The period t osc of small oscillations around the minimum of the potential well 
with V a = 300 and V b = is t osc = 0.84. P$ is the tunneling probability with V 1 = V a -V b = 100. 
Effective chemical potential of the initial condensate: /it e // = 27.15 for N = 10 3 , /i e // = 57.43 
for N = 10 4 . N is the number of 23 Na atoms. Scattering length: a s = 30 A. Inset: tunneling 
probability Pr at t=5 as a function of r. Units as in Fig. 1. 
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